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ABSTRACT 

The anisotropic Heisenberg ferromagnet formalism developed pre- 
viously is examined to include an applied magnetic field for the isotropic 
case in the random phase approximation. Thermodynamic quantities 
such as magnetization, susceptibility and the derivative of magnetization 
with respect to temperature are studied near the Curie point. 

INTRODUCTION 

According to the Weiss theory a ferromagnetic material possesses 
an internal field which is proportional to the magnetization. As the tem- 
perature of a ferromagnet is increased, the magnetization decreases 
until, at a temperature known as the Curie temperature T , the material 
becomes paramagnetic. The process by which this change from ferro- 
magnetism to paramagnetism occurs is referred to as the ferromag- 
netic phase transition. Investigations have shown that electrical, mechan- 
ical, and many thermodynamic properties of a material are altered when 
the material undergoes a phase transition. 

The effect of an external magnetic field is twofold: (1) the magnet- 
ization is increased above its zero field value and (2) the critical point 
disappears. Even though increasing temperature tends to destroy spin 
alinement, the field causes some ordering to be present. Thus instead 
of ah abrupt change at the critical temperature there is a more gradual 
transition. 
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This report is concerned with the investigation of the thermody- 
namic properties of the Heisenberg ferromagnet in an external mag- 
netic field. The quantities studies are the magnetization, suscepti- 
bility and the derivative of magnetization with respect to temperature. 
The last of these can be used to calculate entropy and the magnetocaloric 

effect, which are useful for investigating various recently proposed mag- 

i_2 

netic refrigeration systems. 


Model 


The Heisenberg model is based on a solid where the magnetic elec- 
trons are in states localized about the lattice sites with an exchange in- 
teraction taking place between electron pairs. The model does not take 
into account itinerant electrons and is considered valid for insulating 
ferro magnets such as EuO but not for conductors such as iron and nickel. 
However, it appears that this model gives better results for conductors 

at low temperatures than those calculated from a band theory approach 
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as shown by Argyle, Charap and Pugh. It should be noted, however, 
that the Heisenberg model with its surprising success does not take into 
account the spreading of the electronic energy levels into bands. 

The Heisenberg ferromagnet with spin 1/2 was analyzed by Bogolyubov 
and Tyablikov^ using the techniques of double-time, temperature depend- 
ent Green’s functions. An excellent review of Green’s functions and 
Tyablikov’s application of them to ferromagnetism is given by D. N. 
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Zubarev and references contained therein. 

The Hamiltonian for the Heisenberg model is 


H ■ -^b h o l s i-2 J ij[ s ? s ? + 2 K s j' + s i' s f) 


( 1 ) 


where 




(2) 
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jLtg is the Bohr magneton, g is the Lande g factor, H Q is the applied 

magnetic field which is assumed to be along the z-direction, Sj is the 

spin operator for a spin at site i, and is the exchange interaction 

between spins on sites i and j. The sum is carried over all sites in 

the crystal. The exchange interaction is assumed to be a function only 

of the distance between sites. The self- exchange terms such as or 

J.. are zero. 
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Calculation of Thermodynamic Quantities 

Since the thermodynamic quantities of interest are the magnetization 
and entropy of the system, one is interested in correlation functions of 
the form (s“S^ and hence of the Green’s function ((s^; 

Starting with the Green’s function ^S+j S^Jj one can derive the 
correlation function in the random phase approximation. The results 
are: 



where 


E k = 2<S Z )J(0) 


r 

1 

C 


- + MggH 

J(0)j B 0 


(3) 


(4) 


and 


J(k)=^J gm e ik ' ( 8- m ) (5) 

m 

The sum in equation (3) is over all N lattice vectors in the first 
Brillouin zone. 

The applied magnetic field is assumed sufficient to orient the net 
magnetization along the direction of the field. The field also induces 
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additional long-range order. The long-range order is thus due to the 
internal field and to the applied magnetic field. The thermal average 
of the z component of the spin, <S Z ) , is defined as the magnetization 
per site and is a measure of this long-range order. 

The magnetization, S z is calculated from the relation 

<S Z > = — -1 (6) 

2[1 + 2 cp] 

where 



The sum in equation (7) must be evaluated over all values of k in 
the first Brillouin zone of the appropriate lattice. Except at the very 
low- and high- temperature limits, numerical methods are usually used. 
Such numerical solutions are, however, somewhat difficult. One of the 
purpose of this paper is to show that analytical solutions of equation (7) 
are possible. 

For crystals with cubic symmetry, such as bcc, one can replace 

the sum which appears in equation (7) by an integral. Using the same 

fi- R 

techniques developed by Flax and Raich one can obtain the magnet- 
ization: 

<S Z > = (8) 

2A 

where 



2 + Q 



5 


a _ ^B^o _ H’ 

(ID 

[J(0)t] t 


k„T 


T = B 

(12) 

J(0) 


Q-Ji 

(13) 


T 


and K(k) is a complete elliptic integral of first kind with 



Equations (8) and (9) give <S Z ) a function of temperature, magnetic field 
and magnetization. 

Figure 1 shows a plot of <S Z ) as a function of temperature for sev- 
eral values of H’ for a bcc lattice. When H’ equals zero there are 
no solutions for <S Z > above a critical temperature T , that is, the 
Curie temperature. This is the point where long range order disappears 
and above which the spins are completely disordered. For non-zero H’ 
the ferromagnetic transition occurs over a range of temperatures which 
forms a Curie region, so that the transition is smeared out and results 
in the appearance of a ’’tail. ” 

The cause of the ’’tail” is that the long-range order persists through 
the effect of the applied magnetic field. The higher the field, the greater 
is the broadening of the transition. 

Figure 2 shows a plot of <S Z ) as a function of H’ for several values 
of r for a bcc lattice. As is seen from the figure the magnetization has 
a nonlinear dependence on the magnetic field near the Curie point 
(r = 0. 3588). This nonlinearity in the transition region can be attributed 
to the exchange interaction. This behavior persists as the temperature 
is increased for small magnetic field strengths. For high field and 
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high temperature the magnetization curve becomes linear. This has 
been experimentally verified for many ferromagnetic substances. ^ 


MAGNETIC SUSCEPTIBILITY 
The magnetic susceptibility x is defined as 
,z A 


X 


= jO. = Q 

H’ a 


(15) 


Using equations (8) through (14) together with equation (15) one obtains 
1 


X = 


2a r A’ 


(16) 


where 


A’ 

\P/ 


2K(k,f 

2 2 
+ coth p - I + a 2 JL 
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p = a | (1 + 2x) 
2 2 


(17) 


(18) 


and 


k 2 = ^ 1 -./l - 


2x 


(19) 


2 V VI + 2 X > 

For a given ferromagnetic material the spontaneous magnetization 
can occur only below a critical temperature T c> Well above the Curie 
temperature, such materials are found experimentally to behave para- 
magnetically, and have a susceptibility which follows the Curie - Weiss 
law, namely 


X = 


(T- e ) 


( 20 ) 
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in which C is the Curie constant and 6 is called the paramagnetic 
Curie temperature. 0 is usually somewhat higher than T . 

V 

A plot of 1/x as a function of r from equation (16), is shown in 
figure 3. As required the curve becomes linear at temperature suffi- 
ciently far above the Curie point. The nonlinearity in the neighborhood 
of the Curie point agrees qualitatively with experiment. ** Below the 
Curie temperature, T ferromagnetic materials show a marked in- 
crease in susceptibility. 


CALCULATION OF THERMODYNAMIC QUANTITIES 


For an adiabatic process, where there is no change in entropy 

dr = — — j dH’ (21) 

C H’ \ 9t /H’ 

where C H , is the magnetic specific heat. Since (8 S z /3t)jj! is negative 
for a ferro magnet, an increase in field produces an increase in tem- 
perature; moreover, the increase is expected to be largest near the 
critical point. This is known as the magnetocaloric Effect. 

The entropy, which is a measure of the order of the system, can 
be calculated from 


S = 


nr c 


H 


AH’ 


JO 


dr + 

T JO 


a <s z y 

8r 


dH’ 


'H* 


( 22 ) 


The change of magnetization with respect to the reduced temperature at 
constant applied field can be derived from equations (8) and (9) as follows 



1 

2 




(23) 
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(27) 
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and 


E(k) is a complete elliptic integral of the second kind (28) 

Figure 4 shows a plot of (9 <S z )/9t)jj, as a function of temperature 
for several values of H* for a bcc lattice. These derivatives have not 
previously been obtained analytically. 

CONCLUSIONS 

The thermodynamic properties of the Heisenberg ferromagnet in a 
magnetic field have been investigated theoretically by using a random 
phase approximation. Expressions are derived from which the thermo- 
dynamic parameters can be calculated. The derivative of the magnet- 
ization with reduced temperature is obtained for use in the calculation 
of the magnetocaloric effect and entropy. The analysis predicts that: 

1 . In the presence of a magnetic field the ferromagnetic transition 
is smeared out. This results in the appearance of a tail in the mag- 
netization versus temperature curve. Increasing the magnetic field 
increases the broadening of the transition. 

2. There is a nonlinear dependence of the magnetization upon field 
at the Curie point. When the temperature is raised, the nonlinearity 
remains in the weak field region, but in the strong field region the mag- 
netization becomes linear with field. 

3. Well above the Curie point the susceptibility follows the Curie - 
Weiss law (x"^ proportional to temperature). In the immediate neigh- 
borhood of the Curie point, however the dependence of -1 on tempera- 
ture is nonlinear. These results are in qualitative agreement with ex- 
periment. 
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Figure 1. - Magnetiza- Figure 2. - Magnetization as a function of reduced magnetic field, 

tion as a function of 
temperature. 
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Figure 3. - Reciprocal 
of susceptibility as a 
function of reduced 
temperature. 
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Figure 4. - Partial derivative of mag- 
netization with respect to reduced 
temperature for constant reduced 
magnetic field as a function of re- 
duced temperature. 
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